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ABSTRACT 


Characteristics  of  periodically  stiffened  membrane  are 
studied  by  using  a string  model  loaded  with  equally-spaced  mass 
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Two  theoretical  predictions  of  radiated  sound  power  are 
made  using  statistical  energy  methods.  One  is  based  upon 
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response  of  the  membrane.  Both  agree  well  with  direct 
measurements  of  radiated  sound  power. 
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NOr^NCLATURE 


b Dimension  of  the  cross  section  of  stiffener 

C „ Membrane  wave  speed  for  (m.n)  mode 
m#n  r 

CQ  Acoustic  wave  velocity 

f Frequency 

k^  Wave  vector  in  ^-direction 

k^  Wave  vector  in  ^-direction 

Spacing  between  stiffeners 

Transverse  membrane  dimension 

m Membrane  area  density 

N Number  of  stiffeners  plus  one  (=10) 

n Frequency  modal  density 

s 

p Wall  pressure 

P Acoustic  pressure  radiated  by  membrane 

T*  Membrane  tension 

U Convection  velocity 

Free  stream  velocity 
Z Radiation  impedance 

<J>  (cu)  Wall  pressure  spectral  density 

n Loss  factor 

p Fluid  (air)  density 

tt  Spectral  density  of  radiated  acoustic  power 

a Membrane  area  density 

u Angular  frequency 
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INTRODUCTION 


Much  work  has  been  done  on  the  vibratory  response  of 
thin  plates  or  membranes  to  a turbulent  boundary  layer  pres- 
sure excitation  and  the  subsequent  radiation  of  sound  by 
these  plates  or  membranes . In  the  previous  work,  studies 
have  included  only  the  response  of,  or  radiation  from,  simple 
structures.  Actual  structures  may  be  a complex  system  of 
plates  forming  a skin  over  a rib-like  structure  of  beams. 

In  one  limit,  very  rigid  beams  divide  the  skin  into  a number 
of  adjacent  bays:  the  response  of  such  structures  can  be 

modeled  by  treating  the  response  of  each  bay  separately. 
Generally,  however,  the  beams  are  not  infinitely  rigid,  and 
we  must  consider  the  combined  response  of  ribs  and  skins. 

To  understand  the  response  of  such  systems,  in  the  present 
study  we  investigate  the  response  and  radiation  of  a periodi- 
cally stiffened  membrane. 

The  characteristics  of  the  membrane  and  stiffeners  sys- 
tem are  studied  by  using  a one-dimensional  string  model  car- 
rying equally-spaced  blocks  having  mass  as  well  as  rotary 
inertia.  The  analysis  used  here  is  similar  to  that  used  by 
J.  W.  Miles. * The  natural  frequencies  fall  in  periodically- 
spaced  groups,  the  lower  ends  of  which  correspond,  at  high 
frequency,  with  the  natural  frequencies  of  each  of  the  string 
segments  between  the  blocks  and  with  the  blocks  fixed  as 
rigid  supports . 

The  radiated  acoustic  power  from  the  periodically  stif- 
fened membrane  excited  by  turbulent  boundary  layer  is 
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estimated  using  a modal  analysis  on  an  individual  bay  of 
membrane  between  stiffeners.  The  contributions  of  Lyon  and 
Maidanik/  Leehey,  Davies,  and  many  others,  have  provided 
a broad  fundamental  theoretical  basis  upon  which  this  part 
of  the  study  develops.  The  analysis  is  based  on  an  expan- 
sion of  the  transverse  velocity  field  of  the  membrane  bay  in 
terms  of  orthogonal  characteristic  functions.  It  is  assumed 
throughout  the  analysis  that  the  response  of  the  membrane  at 
any  frequency  is  determined  by  the  response  of  those  modes 
whicn  are  resonantly  excited  at,  or  near  to,  that  frequency. 
This  assumption  is  justified  for  the  cases  of  turbulent 
boundary  layer  excited  structures  when  the  correlation 
lenaths  of  the  exciting  field  are  small  compared  to  the 
panel  dimension.  Then,  expressions  for  the  acoustic  power 
radiated  in  narrow  bands  of  frequencies  can  be  obtained  by 
summing  over  those  modes  that  have  resonant  frequencies  with- 
in the  banc  by  making  use  of  the  statistical  energy  analy- 
sis (SEA)  assumption  of  equal  energy  among  the  modes. 

Corcos's  model  of  the  wall-pressure  correlation  func- 
tion is  used  to  obtain  the  modal  wall-pressure  spectral  den- 
sities. Values  of  the  modal  radiation  coefficients  for 

4 

plates  have  been  calculated  by  Davies.  His  results  are  ap- 
plied  here  to  the  cases  of  membranes. 


II.  CHARACTERISTICS  OF  PERIODICALLY-STIFFENED  MEMBRANE 


The  equation  of  motion  for  a membrane  can  be  written  in 
the  form 

^2  9* 


a2h 


l — 7 f — 7)  h “ f 

zx  az  at 


where  Tmis  membrane  tension,  o the  area  mass  density  of  the 
membrane,  h(x,z,t)  the  normal  displacement  of  the  membrane, 
x-coordinate  is  in  the  longitudinal  direction  of  membrane, 
z-coordinate  is  in  the  lateral  direction  of  membrane,  and 
is  parallel  with  the  axes  of  the  stiffeners. 

Note  that,  if  we  let  h assume  a form  of  h=X(x) Z (z) elwt 
where  X and  Z are  respectively  functions  of  x and  z only 
then  we  can  easily  decompose  the  equation  of  motion  into  two 
uncoupled  one-dimensional  equations: 

afs  ♦ k?  x = o 

3x^  x 


| + k^  Z = 0 


where  kj^  + k^  = k2  = cj2/C2  • cin  being  the 

membrane  wave  speed. 

By  tapering  the  ends  of  the  stiffeners,  and  by  knowing 
that  the  levels  of  vibrations  and  rotations  of  the  membrane 
system  are  small,  we  can  approximate  the  longitudinal  charac- 
teristics of  this  system  by  that  of  a string,  carrying  equal- 
ly-spaced blocks  of  mass  and  rotary  inertia  resting  on  elastic 
foundation.  For  simplicity,  the  elastic  foundation  is  repre- 
sented by  two  springs,  each  of  stiffness  S/2,  for  each  block  as 
shown  in  Figure  2.1.  The  equation  of  motion  of  the  string  is 


f-J  ♦ *!  * - o 

3x 


(2.1) 


where  we  use  y(x,t)  as  the  displacement  of  the  string. 

In  the  formulation  of  the  boundary  value  problem  for  this 
periodical  structure,  we  also  find  it  convenient  to  introduce 
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a dimensionless,  local  coordinate  such  that 

y — ^2 yn ( c ) , ^ = 1,2,.. .,n 


r = 


x-dn 


, d„  < X <■  d„  + i and  d„  = (n-DU.+b) 
n — — n n i 


(2.2a) 

(2.2b) 


where  N is  the  number  of  string  segments  separated  by  stif- 
feners, i.e.,  the  number  of  stiffeners  plus  one,  is  the 
spacing  between  stiffeners.  In  addition,  let 

(2.3) 


A2  = k2i2 


Note  the  A is  also  a dimensionless  eigenvalue  parameter. 
Now,  Eg. (2.1)  becomes 


a2y 


35‘ 


n + A2  y = 0 


n 1,2,.. .,N 


(2.4) 


3y  (0)  <>y  . (1)  ^ i 1 

(~St S1— ) - T^IW11  + Vn<0»j  (2.5a) 


Let  us  study  the  equilibrium  conditions  associated  with 
block  number  (n-1) . We  assume  that  all  angles  are  very  small. 
Equilibrium  of  force  in  vertical  direction  is  thus 

2. 

' J n-1 -i  _A 
35  “35 

where  M is  the  mass  of  one  stiffener,  T = 1^*3,  and 

n ~ 2,3,...,N. 

Assuming  0^  very  small,  the  equilibrium  of  momentum 
with  respect  to  the  center  of  the  block  can  be  written  in  the 
form 


-T(b  sin  0^+b  cos  0^)  + 


/si-1  °rn-l  \ 

(-J-»-l‘1,+T  -M— )(b 

ay„(0) 


co h 0fa-b  sin  0^) 


+T(-b  sin  G^+b  cos  0^)  + 


or 


t si  1 ay„'u;\ 

I- -j- yn  ( 0 ) +T — ^■•)(b  cos  eb+b  sm  %)=lG|i5 


3y  ,(!)  3y.(»)  , . /3V0)  3)rn-lU) 

—Sr sin  6b)  + (—57 sin  °b)  + sin  °b  ( 5 X 57“ 

SSL . 10. 

- sin  0b  (|6  + yn-!(i)  + yn«»)  - fT 


-5- 


where  2b  is  the  dimension  of  the  block.  On  the  left-hand  side 
of  the  above  equation,  the  third  term  is  much  smaller  in  magni 
tude  than  the  first  two  terms,  and,  hence,  is  negligible. 

Also,  remember  that  0b  is  very  small  Now,  the  equation  goes 
over  to 


-20b  + 


Syn-Ia)  3yn(0) 


TT 


TT 


. Sb 
)-  T eb 


i 

tF  °b 


n = 2,3, ... ,N 


Since  the  forcing  functions  of  the  above  differential  equa- 


tions, 3yn_i ( 1) /3 C and  9yn (o)/3C# have  harmonic  time  dependence 
“ 2 

of  0),  0b  can  be  replaced  by  -w  Gb.  Also,  note  that 
yn (°)  ” yn-l ^ * I^^b  ; n = 2»3'***N*  Thus,  we  have 


+ V°>  ' 


2b  . 3V°> 

T 2 ,b  l SC  + 4c 
(2-^-#)  *x 

1 n = 2,3,. ..,N 

On  the  ends  of  the  string,  we  have 

y1 (0)  * 0 

yNd)  * o 

The  most  general  solution  to  Eq.(2.4)  is 
yn(S)  = An  sin  A£  + cos  A £ ; n = 1,2,...,I. 


) 


(2.5b) 


(2.5c) 

(2.5d) 


(2.6) 


Now,  Eq. (2.5a)  and  Eq.(2.5b)  give,  respectively 

ct[B  +A  , sin  A+b  . cos  A]  = A_  - A^  . cos  A 4-  B . sin  A 

n n-i  n-i  n n-1  n-1 


3 [B  -A  . sin  A-B  . cos  A]  = A + A . cos  A - B , sin 
n n-i  n-l  n n-1  n-1 


where 


(2.7a) 


(2.7b) 


a = 


3 = 


(S-Mu>  ) iL1 
2TX 


bA 

n = 2|3|«t«fN 


Eq.  (2.7a)  subtracted  from  (2.7b)  gives 


- (a-$)B  + [- (a+0)  cos  A+2  sinX  ]B 


From  Eqs. <2.10a)  and  (2.10b),  we  have 

(a-g)  sin  A sin (N-2) 0+2 (cos  A+a  sin  A) (cos  A+g  sin  A) sin (N-l) 0=0 

(2.11) 

From  Eq.(2.10),  we  see  Z = cos  0,  where  0 is  real  when 
| Z | is  equal  or  less  than  one,  and  0 is  imaginary  when  |z| 
is  greater  than  one.  It  is  easy  to  see  that  when  0 is  real, 
we  will  be  able  to  solve  for  eigenvalues  from  Eq.(2.11).  If 
0 is  imaginary,  say,  0 = ix,  with  X real,  Eq. (2.11)  becomes 
(a-g)sin  A sinh (N-2) X+2 (cos  A+a  sin  A) (cos  A+g  sin  A)sinh(N-l)X 

= 0 (2.11b) 

A little  calculation  shows  that,  for  large  value  of 

A, (a-B)  sin  A is  less  than  2 (cos  A+a  sin  A)  (cos  A + g sin  A) . 

Besides,  we  know  sinh (N-l) X is  larger  than  sinh(N-2)X.  Thus, 

it  is  impossible  to  find  a solution  of  A which  makes  ] Z | >1 

and,  at  the  same  time,  satisfies  Eq. (2.11b).  So,  we  see  the 

resonant  frequencies  are  permitted  only  within  bands  where 

|z|<_l  or  0 is  real,  as  shown  in  Figure  (2.2).  Note  that  at 

high  frequency  the  asymptotic  approximation  of  Z is 

2 2 

Z a,  c A sin  A,  where  C = IM/ (M  - I/b  ) loJt^b  . A computer 
analysis  shows  the  resonant  frequencies  as  listed  in  Table 
2.1.  Notice  that  in  each  frequency  band  there  are  ten  reso- 
nant modes.  We  applied  Lagrangian  energy  method  to  the  system 
with  a simplification  that  the  strings  are  massless.  We  found 
two  natural  frequency  bands  comparable  to  the  first  two  bands 
listed  in  Table  2.1,  and  that  the  first  band  involves  mainly 
the  heaving  and  the  second  the  rotating  of  the  blocks.  This 
is  a rather  comforting  check  of  the  analysis  we  pursued.  The 
modal  shape  of  higher  frequency  bands  are  characterized  by 
vibrating  and  rotating  blocks  linked  b',  string  segments  which 
assume  the  modal  shape  of  a string  of  length  whose 
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fundamental  frequency  coincides  the  lower  frequency  limit 
of  each  band.  Also  notice  that  as  frequency  increases  the 
bandwidth  of  each  group  of  modes  decreases. 

Table  2.1 

String  Modal  Resonant  Frequencies  (Hz) , Determined  by 

Eq. (2.11) 


Band  1 
Band  2 
Band  3 
Band  4 
Band  5 


222  225  231  241  253  267  282  295  306  313 

840  846  862  886  917  952  989  1026  1057  1079 

1525  1547  1577  1610  1645  1677  1704  1725  1738  1740 

2999  3006  3016  3030  3044  3058  3071  3081  3087  3090 

4485  4489  4496  4504  4513  4522  4531  4537  4541  4543 


Now  we  return  to  the  two-dimensional  membrane.  The 

modes  of  the  membrane  are  those  of  the  above-analyzed  string 

coupled  with  those  of  a string  of  length  the  width  of 

the  membrane,  i.e.,  (k  ,k3)  where  k1#k3  are  respectively 

wave  vectors  of  the  abovemen tioned  string  systems.  This  is 

2 2 2 2 

best  shown  in  Figure  (2.3).  Note  that  k * k^  + * w a/T, 

and  that  f = 2-rrkC^,  where  is  independent  of  fre- 

quency, i.e.,  f varies  linearly  with  k.  Since  in  k-space 
the  distance  from  origin  to  mode  (k^,k3)  is  just  k,  we  can 
easily  transform  modes  from  ic-space  to  frequency  domain. 
First.,  we  use  a circle  centering  at  the  origin  to  pass  the 
mode  concerned  in  two-dimensional  ic-space  and  cross  the 

k-scale  at  k which,  aft^r  multiplied  by  2irC  , is  just  the 
p • m 

resonant  frequency  of  that  mode  (Figure  2.3). 

III.  THE  RADIATED  POWER 

A high  frequency,  the  stiffeners  on  the  membrane  act 
as  rigid  barriers.  This  explains  why  the  ic-space  picture  of 
the  membrane-stiffeners  system  looks  much  like  that  of  an 
individual  membrane  bay. 
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In  evaluating  the  radiated  power  of  this  system,  we 
made  two  assumptions: 

1.  The  radiation  contributions  from  different 
membrane  units  are  statistically  uncorrelated. 

2.  Each  membrane  unit  radiates  as  if  it  were 
mounted  in  ar.  infinite  rigid  baffle. 

These  two  assumptions  need  discussion.  The  first  one 

is  justified  because  the  exciting  wall-pressure  field  is  not 

correlated  between  bays  (cf . Section  V) . 

Now  let  us  consider  the  following  case:  A membrane  is 

excited  on  one  side  by  a pressure  field,  P,  and  is  mounted 

in  an  array  of  membranes  of  the  same  size.  The  acoustic 

power  radiated  into  the  other  side  of  the  membrane  includes 

2 

the  direct  contribution,  <p  >,  from  the  excited  membrane  plus 
the  back  reaction,  <p*',  from  the  rest  of  the  membranes  ex- 
cited by  p.  Let  °rad  be  the  radiation  coefficient,  and  Z the 
total  damping,  i.e.,  both  mechanical  and  radiation,  of  the 
membrane.  If  v is  the  velocity  field  of  the  membrane,  then 
<v  > "v  <P  >/|z|  . Also,  we  know  <p  >/pc  *v»  ir  ^ % °rad  <v  > 

where  is  the  acoustic  power  radiated  directly  from  the 

2 2 2 

membrane.  Thus,  <p  > ^ <P  > Pcc>rad/|z|  . in  general, 

2 2 

<p  > <<  <P  >.  The  same  argument  leads  us  to  the  conclusion 
2 2 

that  <p  > <<  <p  >.  And  this  justifies  assumption  2 stated 
above . 

After  the  appropriateness  of  these  two  assumptions,  we 
can  consider  each  bay  as  an  independent  radiator,  and  fur- 
ther approximate  each  radiator  as  a rectangular  membrane 
unit  of  length  in  the  direction  of  mean  flow  and  of  width 
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l 2 inserted  in  a rigid  baffle.  The  acoustic  radiation  into 
the  side  away  from  the  flow  is  considered  here.  The  analy- 
sis given  in  this  and  the  following  two  sections  is  a repeat 
of  Davies'  work  on  rectangular  plates,  and  is  beirg  applied 
to  each  bay  of  our  membrane  and  stiffeners  sysJ  . high 
frequency . 

We  choose  Cartesian  Coordinates  with  x^,x^  in  the  direc- 
tions of  and  l 3 respectively,  and  x2  normal  to  the  mem- 
brane, but  away  from  the  flow  side.  The  equation  of  motion 
for  membrane  displacement  y(x,t)  driven  by  a wall  pressure, 
p (x , t ) , is 


TV2y  + mg  |£  + m — £ = -p(x,t)  + P(x,x2=0,t)  , 

3 1 


(3.1) 


where  T is  the  tension  of  the  membrane,  m is  the  area  den- 


sity of  the  membrane,  8 = n |w|  accounts  for  mechanical  dam- 

s 


ping  of  the  membrane  through  a linear  dependence  on  velocity. 


's 


is  the  loss  factor,  and  P(x,x2,t)  is  the  acoustic  pres- 


sure generated  by  the  motion  of  the  membrane. 

At  the  interface  x2  = 0 , the  boundary  condition  is 


SP 

3x. 


3v 


= -P  TT 


3t 


(3.2) 


where  p is  the  fluid  density,  and  v(x,t)  = 3y/3t.  The  dis- 
placement y vanishes  at  the  edge  of  the  bay. 

The  transverse  velocity  field  of  the  membrane  can  be 
expanded  in  terms  of  orthogonal  characteristic  functions 
(the  in  vacuo  normal  modes) . The  frequency  Fourier  trans- 
form of  membrane  velocity  is  written 

v (x , cj ) = /*  v(x,t)eiujt  dt  = l 


v ( O ) * (X) 

_ _ , mn  mn 
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where  ij>  (x)  = ( 2//A ) sin  k x.  sin  k x0.  Here  k = mir/SL.  , 
inn  in  x n j in  j. 

kft  = and  A = is  the  area  of  the  membrane  unit. 

The  modal  equation  obtained  from  Eq.(3.1)  has  the  form 


2 2 2 

(-T  k + iw  mn  sgn  w-to  m)  v ( cu ) = 

ion  s mn 

+ iwp  (w)  - iwP  (w)  , 

rmn  ' mn  ' 

2 2 2 

where  = km  + kR,  and  and  Pinn  are  defined  by 

*nn(w)  = / P(x/W)  ip  (x)  dx, 

A 

Pmn(uj)  = / P(x'w)  ^mn  (^J  dx  * 

A 


(3.3) 


The  wavenumber-frequency  Fourier  transforms  of  acoustic  pres- 
sure and  membrane  velocity  are  connected  by  using  Eq.(3.2), 
together  with  the  wave  equation  for  the  acoustic  medium. 

Thus , 

P(k;::2,w)  = 2 (it  ,w)  V (it, a;)  exp [ix2  (k2-k2) 1/2]  (3.4) 

where  it  is  the  wavenumber  vector  on  the  membrane,  Z(it,w)  is 
the  radiation  impedance  defined  by 


Z(it,w)  = p C (1  - k2/k2)“1/2  , 

00  o 

where  k = |k|,  and  kQ  = w/Co'  Co  keing  t^le  wave  velocity  in 
the  acoustic  medium. 


Now  P (■>)  can  be  written  in  the  form 
mn 


P (uj)  - 
mn  ' 


(2tt)  qp 


I%p(w)  r r z<£,u>)  snn  (it)  s*ti  (it)  dit  , 


—00  —00 


qp 


mn 


(3.5) 


where  S (it)  = f (x)e^k’xdx 
mn  ^ rmn 


and  the  asterisk  denotes  a complex  conjugate. 

Ve  consider  here  a system  where  the  response  is  prima- 
rily due  to  resonant  modes  and  are  thus  interested  in  the 

value  of  P only  at  the  resonant  frequency  of  the  (m,n) 
mn 
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modes.  If  the  resonant  frequencies  are  widely  separated  in 

frequency#  (nt/n)  is  the  only  effective  terra  in  the  summation. 

Pran((°)  t*lus  can  approximated  by  one  terra 

P_(m)  = P C v (ui)  (a  ix__) 

ran  o o ran  ran  ran 

With 


a + iy 
ran  *mn 


jHr  / z<*'*)lsran(*)|2  ^ 


poCo 


(3.6) 


With  this  approximation,  the  modal  velocity  is  given  in 
terms  of  the  modal  wall  pressure  p^  by  the  equation 

v (w)  = Y (to)  p (w)  , 

tun  ran  ran 

From  Eq. (3.3),  is  obtained. 


(3.7) 


Ymn(w)  “ iaj ( ^“Tkinn”aj  m"wpoCoxran)  + i(ui  ms  sgn  w+upoCo0inn;  ) 


(3.8) 

It  is  shown  in  Section  II  that  the  virtual  mass  term  invol- 
ving  is.  in  general,  negligible  for  acoustically  slow 
modes  if  the  acoustic  iridium  is  air. 

The  wavenumber-frequency  transform  of  acouptic  pressure 
can  now  be  written  in  the  form 

P(k,x2,oo)  =Z(£,uj)  l Ymntw)Pnin<w)  ^(^J  exP  Ux2  (k^-k2)1/z] 

w,n 

And  P(x,x2,oj)  can  be  obtained  from  the  inverse  transform. 

Thus,  the  spectral  density  of  the  power-radiated  tt  (w)  , is 
written 


Tf(w)  = / E[P  (x,o,uj)  v*(x,w)J  dx  , (3.9) 

A 

where  it  is  understood  that  the  real  part  of  the  integral  is 
taken  as  the  acoustic  radiation. 

Substitute  the  expressions  of  panel  velocity  and  acous- 
tic pressure  into  Eq.(3.9)  and  we  have  the  spectral  density 
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writ  ten  in  the  form 
1 


tt(w)  - 


4 


l l r Z(it/W)smn(^)s*n(}t)d}c 


( 2-TT > m,n  q,p  -« 


qp 


X Y,nn(“)Yqp<“>  [I  ♦<£.“>Smn<*>S5p(S)dS  , 

where  <$>(&/“)  is  the  wavenumber- frequency  transform  of  the 

wall-pressure  cross-correlation. 

7 

Using  Dyer's  criteria,  we  can  reduce  the  summations  of 
Eq.(l.lO)  to  a single  infinite  summation.  Essentially,  it 
is  assumed  that  4>(Jc,w)  is  constant  in  Jc.  The  cross  terms 
in  Eq.(3.10)  thus  vanish  because  of  the  orthogonality  of  the 
modes . 


»oCo  Z 
m,n 


where 


* 77TT  I |Smn<r<>|2  di 

( 2t ) 


(3.11) 


(3.12) 


:'.s  the  modal  wall-pressure  spectral  density. 

Again,  since  we  are  treating  a highly  resonant  system, 
the  radiated  power  in  a narrow  frequency  band  is  dominated 
by  radiation  from  those  membrane  modes  that  have  resonant 
frequencies  lying  in  the  band.  Eq.(3.11)  can  then  be  inte- 
grated over  the  narrow  band  Aw  centered  on  wq  to  give 


tt  (w  ) Aw  = p C y a (w„)  <p  (w  ) /° 

' r\ ' ^ r\  r\  mn  ' o Tmn  * 

m,n 


w +Aw/2 


Aw, 

uo”  '2 


dw  , 


wheie  the  summation  now  includes  only  those  modes  resonant 

in  the  Aw  band.  The  in  egrai  is  *vxaluated  by  writing 

w = w (1+Y)  - w (1+Y),  where  Y <<  1 for  frequencies  in  the 
mn  o 

Aw  band.  Also,  let  us  define  = ns  + (p0cG/U)0mp)  °mn*  a 

loss  factor  including  both  mechanical  and  acoustic  losses. 


& 3>  s*"**«3F*5r?-‘ v ' 


:3ssn 
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Fur  ther.  we  assume  n «<<  Aw/2u)  <<  1.  Thus, 

mn  o 

m+Am/2  ~ Aw/2oj  ~ 

Lw2|V”»("H  d“*2£  °“o|Y-l<,lf 

And  now. 


IT 


2Vom 


71  = PcCo 


T 


( I 


(3.13) 


2wm  m,n  'ran 

In  the  present  experiment,  it  was  only  possible  to  mea- 
sure total  loss  factors  in  1/3-octave  frequency  bands.  Thus, 

it  is  reasonable  to  assume  a constant  n „ for  the  frequency 

mn 

band.  Also  note  that  4>  /n  „ is  proportional  to  the  modal 

mn  mn 

energy.  An  assumption  of  statistical  energy  analysis  (SEA) 
is  that  the  power  flow  from  the  membrane  to  the  acoustic 
medium  can  be  estimated  by  treating  the  modal  energies  of 
all  modes  in  the  frequency  band  as  equal.  Then  can 

be  taken  outside  the  summation  in  Eq. (3.13)  and  the  expres- 
sion rewritten  in  the  form 


tt(oj)  = p C (7t/2nura  ) n > 

O v O 


mn  ’mn 

where  <•••>  represents  an  average  value  over  reson-m  modes 

in  the  Am  band  and  n is  the  frequency  mode  density. 

s 

IV.  MODAL  RADIATION  COEFFICIENTS 

Eq. (3.6)  can  be  written  in  the  form 
16  k 


(3.14) 


a + i y 
mn  Amn 


r f 
<2tt)  A o o 


I.(k.)I.(k.) 

■ 2 2 TT77  ^^1^3  * 

(k^-k‘-k“)v<£  1 i 


(4.1/ 


where 


2 2(l*(-l)m  cos  k^J 


h = Z2 

*m 


5 — JT 

u-kf/kfV 

I m 


(4.2) 


and  I ^ has  a similar  form.  The  membrane  studied  in  the 
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present  report  has  a wave  speed  of  about  1/3  the  acoustic 

speed  in  the  air.  Thus,  the  functions  here  are  evaluated 

only  for  acoustically  slow  modes;  in  fact,  the  restriction 

k^  > 2kQ  is  imposed.  Therefore,  the  integral  for  always 

includes  the  point  k.  « k , k,  = k . Then,  both  I,  and  I, 

i m j n i j 

can  be  approximated  by  delta  functions  (see  Reference  6) . 


s (kx-V 


The  integration  now  yields  ^ ~ kc/kmn*  11116  v>irtual 
term  in  Eq. (3.8)  is 

p C _ p . _ 

u)2m(-~-^  Xmn)  = w2m(-2  < ui  m (0.05)  . 

aim  ran  m K 

mn 

where  density  of  air,  p , is  1.29  Kg/m3,  surface  density  of 

2 

the  membrane,  m,  is  0.292  Kg/m  , and  k^  is  given  its  smallest 
value  in  the  experiment  reported  here,  80  1/n  (at  1500  Hz). 

The  virtual  mass  term  can  thus  be  neglected. 

Approximations  to  the  modal  radiation  coefficients 
for  various  frequency  regimes  have  been  given  by  Maidanik 
and  Davies.  The  results  are  applied  here  to  our  cases. 


First,  consider  cases  where  k > 2k  . 

mo 

(a)  n = 1 or  2 

Approximate  results  are: 

. k £ sin(k  L/2.4) 

1 32tt  o 3 r-i  , 0 3 

°ml  " 7TTZ  ~ 7TT  ^ TT77T71 


- (-1) 


_ sin  k i- 
rn o i 


k l&2/(2.4) 2 + i?]1/2 


k U:/(2.4)“  + £7 


1 k sin(k  JL./1.7)  _ sm  k 

■*  Tw7  ~ 11  1 ' ~W1.>  ‘-1* 

sin  k nl/a.7)2  + a?] 1/2 

+ (-Dm r°-  ^ ) (4.3) 

k0U3/(1.7r  + 

(b)  kn  < kQ  (m,n>2) , edge  modes 

ran  ‘ -~7  i1  - (-1)m  Wko  - kn)1/2li  <4-4’ 

<2ll)  km 

(c)  kQ  < kn  < 2kQ  (m,n>2) , transition  region  between 
dge  and  corner  inodes 


1 2rr 


(2p) 2 A k 


V f.,1  + VS  . 2ko/kn 

:2k  1 “ ko/kn  (1  - k2/k 

m n ' o 


(d)  kn  > 2kQ  (m#n>2) , corner 


ft  - k>n> 
modes 


sin  k £, 


sin  k JL 


+ (-D 


m+n  Sln  ko(fcl  - ‘31 

V*1  + 


Similar  expressions  exist  for  km  > 2kQ  cases. 

The  average  value  of  Eq.(4.6)  can  be  taken  as  represen 

tative  of  the  radiation  coefficient  for  small  values  of 

k £,  and  k and  is  written 

o 1 0 3 

k2 

- 1 32n  *o  ,4 

°mn  “ (7  . 2 A .2$  ' 

(2tt)  kmKn 

Inspection  of  Figure  (4.1)  shows  that  this  typical  value  is 
valid  approximately  for  all  modes  such  that  kQJl^  ko^3  < 
3tt.  The  total  average  value  of  the  radiation  coefficients 


-17- 


of  the  (corner)  modes  resonant  in  the  Aw  band  is  then  ob- 
tained, for  k , k £,0  < 3tt,  as 
o l o 3 


<o  > * -r—  I2  2kpll  o k do  - -|2_  ^ (1,  + l ) 
mn  pk“  J „ *■  mn  P „3A  k3  1 3' 


(4.8) 


For  either  or  both  of  k L or  k £_  > 3tt,  some  modes  con- 

o l o 3 — 

tribute  as  edge  modes.  The  average  value  is  then  obtained  as 

to  K a 

<c  > B o + 2-  o (t  + £ ) . (4.9) 

mn  ttA  ^3  uA  ^3  1 3 

P P 


V.  MODAL  WALL-PRESSURE  SPECTRUM 

Assuming  Corcos's  model  for  the  wall-pressure  correla- 
tion function. 


wr. 


wr. 


, I 1 I I ****  T I la?ri  . 

= *<“>  ex?  (-  aru—  - aTtn o— 0 - 

1 C 3C  c 


where  U is  the  convection  velocity,  and  has  a value  of 
c 

about  0.71)^,  and  <f»(w)  is  the  frequency  spectral  density, 

Q 

Blake  showed  that  = 8 and  = 1.1,  for  the  wind  tunnel 
where  the  experiments  reported  here  were  done. 

The  correlation  length  in  the  x1-direction  is  a^Uc/w, 
according  to  the  above  Corcos's  model.  For  frequency  range 
higher  than  1 kHz,  and  U =48  m/sec,  a.U  /w  is  smaller  than 
Jl^,  the  longitudinal  dimension  of  one  membrane  unit.  This 
justifies  assumption  2 made  in  Section  III.  The  radiation 
contributions  from  different  membrane  units  are  statistically 
uncorrelated. 

4 

It  was  shown  by  Davies  that  Eq. (3.12)  can  be  rewritten 
in  the  form. 


<V“/v  + pi  (kmwuc}  + k;  + 35 


(5.3) 


where  8^  = co/ot^Uc  and  8^  - w/a^Uc.  Thus,  the  total  average 
value  of  the  modal  wall-pressure  spectral  density  over  the 
resonant  modes  in  a A:,;  band  is  then  in  the  form 


IT 

< (j>  > = f (u)  k dO 

mn  uk  J mn  ' p 
p o 


4 ot  ct  U 

= ♦(«)  — -|-i.  jj!  r/2  ( [1+a2  (1-Y  cos  0)  2]“1  + [1+a2  (1+y  cos  0)2]"1) 


x (1  + a2  y2  sin2  0)"1  dO 


(5.4) 


where  7 = U k /oj  = u /C  . 

c p c p 

Since  y < 0.1,  the  integral  in  Eg. (5.4)  can  be  approxi- 


mated as 


<Ln>  = <P  (w) 
ran 


4“3Uc 


(5.5) 


VI . EXPERIMENTAL  RESULTS 


Experimental  Apparatus 


A 28-inch  diameter  mylar  (polyethylene  terephlatate) 
bass  drum  head  was  uniformly  tensioned  on  a bass  drum  frame. 
The  uniformity  of  the  tension  was  determined  by  exciting 
the  circular  mylar  membrane  acoustically  from  underneath 
with  pure  tones  at  its  resonant  frequencies.  The  resultant 
Chladni  patterns  were  observed  using  a light  dusting  of  fine 
sand  to  identify  nodes. 

With  uniform  tension  achieved,  a backing  plate  of  3/4- 
inch  plywood  was  cemented  to  the  lower  surface  of  the  mylar 
membrane.  This  backing  plate  contained  a plexiglass-edged. 
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rectangular-shaped  cutout  which  served  to  define  the  desired 
membrane.  Again,  acoustic  excitation  and  sand  were  used  to 
observe  the  Chladni  patterns  of  the  rectangular  membrane. 

It  was  from  these  observations  that  modal  resonant  frequen- 
cy response  could  be  correlated  to  specific  modes. 

Then  we  cemented  nine  a]  minum  square  bars  periodically 
spaced  onto  the  lower  surface  of  the  rectangular  membrane. 

To  avoid  extreme  local  stresses  due  to  the  square  shoulders 
of  the  stiffeners,  the  bars  are  tapered  at  the  end.  Also  to 
allow  heaving  freedom  of  the  bars  on  the  membrane,  there  are 
small  gaps  between  the  bar-ends  and  the  boundaries  (plexi- 
glass) of  the  membrane.  The  effects  of  these  small  gaps, 
e.g.,  scattering  of  waves  in  the  membrane,  are  not  clear  at 
the  moment.  The  membrane  and  stiffeners  system  had  the 
following  characteristics  (Figure  6.1): 

1.  Membrane  dimensions:  L1  = 54.5  cm,  L3  = 11.5  cm 

2.  Thickness:  8.5  mils 

2 

3.  Surface  density:  o = 0.292  Kg/m 

4.  Spacing  of  stiffeners:  4.59  cm 

5.  Stiffener  dimensions:  10  cm  x .933  cm  x .933  cm 

Before  the  rectangular  membrane  was  stiffened,  for  an 

(m,n)  mode,  its  resonant  frequency  f and  the  membrane 

in  f n 

9 

wave  speed  C „ are  related  by  Morse. 
c m,n  J 


m,n 


\ * C 


m,n 


(r-)2  ♦ (r~) 


where  C"  „ = T/o,  . . , and  a,  , , = o + added  mass.  At  high 
m,n  total  total 

frequencies,  the  modal  added  mass  of  a resonantly  responding 
mode  becomes  vanishingly  small,  hence  the  value  of  C 


m,n 
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approaches  asymptotically  the  in  vacuo  membrane  wave  speed, 

, which  determines  the  membrane  tension,  T.  Before  stif- 
fened by  the  aluminum  bars,  the  of  this  membrane  was 
137  m/sec  (Figure  6.2). 

Precaution  was  taken  to  eliminate  most  vibration  of  the 
backing  plywood  plate  by  attaching,  with  screws  and  cement, 
additional  damping  sheet  lead  1/4-inch  thick  onto  the  ply- 
wood. The  backed  membrane,  still  held  by  the  tensioning 
frame,  was  flush-mounted  in  the  wind  tunnel  test  section, 
forming  the  lower  panel  of  the  test  section,  with  the  longi- 
tudinal direction  in  the  flow  direction.  The  details  of 
membrane  mounting  and  test  section  are  shown  in  Figures 
(6.3)  and  (6.4). 

Total  Damping  of  the  Membrane  and  Stiffeners  System 

These  measurements  were  performed  in  the  following  way. 
We  first  imposed  a D.C.  voltage  to  a shaker  which  would 
push  upward  against  the  membrane.  Secondly,  we  fed  the 
shaker,  in  addition,  an  A.C.  signal  of  1/3  octave  white 
noise,  but  with  the  shaker  tip  in  good  contact  with  the  mem- 
brane. Finally,  we  shut  off  the  shaker  so  its  tip  was  clear 
and  away  from  the  membrane  and  a fiberoptical  displacement 
gauge  flowed  by  1/3  octave  filter,  preamplifier  ard  graph 
level  recorder  would  observe  and  record  the  decay  rate.  The 
reverberation  time,  TR,  is  the  time  needed  for  the  level  of 
a transient  signal  to  decay  60  dB.  Its  relation  with  the 
loss  factor  is  TR  = 2.2/fQnt  where  f is  center  frequency 
of  the  1/3  octave  band. 

The  loss  factor  for  the  system  is  plotted  in  Figure  6.5. 
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Displacement  Measurements  of  Membrane 

The  turbulent  boundary  layer  pressure  fluctuations  were 
created  by  the  mean  flow  in  the  wind  tunnel  at  speeds  of 
32  m/sec  and  48  m/sec.  Again  we  used  the  fiberoptical  dis- 
placement gauge  and  graph  level  recorder  to  measure  and 
record  the  vibration  levels  at  many  different  positions  of 
the  membrane/  in  1/3  octave  bands.  Finally,  we  averaged  the 
data  for  each  of  the  1/3  octave  bands.  These  values  are 
shown  in  Figure  6.6. 

Damping  of  Room 

Standard  techniques  were  employed  in  measuring  the 
reverberation  time,  in  1/3  octave  bands,  of  the  reverberant 
chamber  which  housed  the  test  section  of  the  wind  tunnel. 
Experimental  results  are  presented  in  Figure  6.5. 

Sound  Pressure  Level  Measurements 

3 

The  volume  of  the  reverberant  chamber  is  893  ft  . 

Again,  standard  techniques  were  used  in  these  measurements. 
With  the  wind  tunnel  running  at  32  m/sec  and  at  48  m/sec, 
three  1-inch  B & K condenser  microphones  were  used  to  mea- 
sure simultaneously  the  sound  levels  at  different  positions, 
and  these  levels  are  averaged  for  each  1/3  octave  band. 

There  is  a minimum  of  3 dB  difference  between  the  sound 
pressure  level  measured  with  the  membrane  system  in  place 
and  the  background  noise  level  measured  with  a sand-loaded 
plywood  box  plugged  onto  the  membrane  system. 

The  measured  and  averaged  sound  pressure  levels  are 


shown  in  Figure  6.7. 
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VII.  ESTIMATES  ON  RADIATED  SOUND  POWER  LEVEL 
(a)  Based  on  Displacement  Measurements 

In  a 1/3  octave  band,  whose  lower  and  upper  limit- 
ing frequencies  are  f^  and  f2  respectively,  the  radiated 
sound  power  in  watts  is 

» ■ {2  Rrad  v<f>  df 
X1 


where  R . 3 p C Ao  , is  the  radiation  resistance 

rad  o o rad 

A is  the  total  vibrating  area  of  the  membrane 

orad  is  the  radiation  coefficient 
2 

V (f ) = (2Trf)  D (f ) is  the  spectral  density  of  velocity 
and  D (f ) is  the  displacement  spectral  density 
We  measured  D(f),  and  from  previous  analysis  we  computed 
orad.  From  these  data  we  obtained  the  radiated  sound  power 
level . 


PWL  = 120  + 10  log1Q  it  , dB  re  10”12  watt 
The  results  are  shown  in  Figure  7.1. 

(b'  Based  on  Sound  Pressure  Level  Measurements 

For  a highly  reverberant  room,  the  sound  power 
level  is  given  by  the  following  expression: 


PWL  = SPL  + 10  log1Q  V + 10  log1Q  D - 47.3 

3 

where  V = 893  ft  is  the  volume  of  the  room 
D = 60/T_  is  the  decay  rate  in  dB/sec 
SPL  is  the  average  sound  pressure  level  in  the  room 
The  radiated  sound  power  levels  calculated  by  this 
method  are  presented  in  Figure  7.1. 
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^c)  Based  on  Wall  Pressure  Data 


Wall  pressure  spectral  density,  <J>(f ) , is  readily 

8 

available  from  Blake's  report.  These  data  are  shown  in 
Figure  7.2.  As  shown  in  that  figure,  we  approximated  <p  (f ) 
of  different  frequency  regions  and  of  different  wind  tunnel 
speeds  by  several  straight  lines.  The  linear  functions  for 
these  lines  are  listed  in  Table  7.1. 

By  inserting  4>(f)  into  Eq.  (3.14)  and  integrating  over  a 
1/3  octave  band,  we  obtained  the  radiated  sound  power  levels 


as  shown  in  Figure  7.1. 


Table  7.1 


Approximations  for  <j>(f),  the  Wall  Pressure  Spectral  Density 


U (-^-) 
« sec 


f < 2239 


f > 2239 


32  4>(f)  = .88  8xl0“3-.229xl0**6f  <J>  (f ) = . 538xl0~3- . 712xl(T7f 

48  <p  (f  ) = .337xl0_2-.105xl0_5f  <j>  (f ) = . 104xl0-2 


Note:  The  unit  for  <Mf)  is  Kg2/m2sec3. 

*2239Hz  is  the  upper  limit  frequency  of  the  33rd  1/3 
octave  band  centering  at  2000  Hz. 


Vlll.  DISCUSSION 

We  sampled  the  vibratory  response  at  many  different 
positions  on  the  periodically  stiffened  membrane.  Figure 
8.1  shows  some  typical  examples.  The  samples  show  sharp  and 
clear  frequency  band-pass  characteristics , as  the  analysis 
in  Section  II  pointed  out. 

We  used  three  methods  to  estimate  the  radiated  sound 
power  levels  at  two  different  wind  tunnel  speeds,  32  m/sec 
and  48  m/sec.  The  first  one  was  based  on  displacement  mea- 
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surements,  the  second  one  on  sound  pressure  level  measure- 
ments/ and  the  third  one  on  wall  pressure  data.  The  experi- 
mentation showed  good  agreement  among  these  three  estimates. 

In  Figure  8.1/  the  first  peak  at  about  250  Hz  is  due  to 
the  "heaving"  modes  or  band  #1,  as  predicted  in  Table  2.1. 

The  second  peak  near  450  Hz  is  due  to  the  "pitching"  modes 
not  covered  by  the  analysis  in  Chapter  II.  This  is  easily 
understood  by  modeling  our  system  as  rectangular  bars  each 
resting  on  two  springs  at  the  ends  of  the  bar,  as  shown  in 
Figure  8.2.  The  ratio  between  the  pitching  and  heaving  reso- 
nant frequencies  is  about  /3,  and  that  is  the  ratio  between 
the  second  and  the  first  peaks,  as  discussed  above. 

We  make  the  following  study  to  see  the  reason  why  band 
#2,  from  840  Hz  to  1079  Hz,  listed  in  Table  2.1,  does  not  ap- 
pear in  the  displacement  spectrums.  First,  remember  that  the 
fundamental  resonant  frequency  for  a x membrane  unit  is 
1525  Hz,  i.e.,  the  lower  frequency  limit  of  band  #3  (see  Fig- 
ure 2.3).  Note  that  band  #1  and  band  #2  are  well  below  that 
frequency.  This  means,  in  the  string  and  block  model,  that 
the  strings  between  the  blocks  are  nearly  straight  segments. 

For  c-uch  a system,  the  cutoff  frequency,  f , is  associated  with 
A = l 1 because,  for  A < £, , there  will  be  cancellation  of  ex- 
citations.  Secondly,  note  that  heaving,  h^  « 0^  - 0^_^  and 
rolling,  A^  11  0j_  * Qi_i  (dee  Figure  8.3).  Also  note  that,  for 

high  frequencies,  the  root  mean  squares  of  heaving  and  rolling, 

~~2' " — “•> — y — 2 y 

h^(oo)  and  t'^(oj)  are  both  proportional  to  |0^(oj)|  + j©f_^(co)|  , 

since  th.5  cross-correlation  of  and  0^_^,  0^  (w)  (oi)  = 0, 

which  is  true  at  least  for  A < So  we  may  conclude  that 


ubmImMium 
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heaving  and  rolling  modes  both  have  the  same  cutoff  frequency. 

Next/  we  know  that,  for  a frozen  convecting  turbulence, 
k = 2tt/A  = 27Tf/U  where  U = . 7U  . The  larger  value  of  tl  is 
48  m/sec  in  our  experiments  and  = 4.6  cm.  Thus,  the  cutoff 
frequency  is  fQ  = 0.7  - 700  Hz.  Finally,  it  is  quite 

clear  now  that  band  #2,  beyond  the  cutoff  frequency,  will  not 
be  excited  by  the  turbulent  boundary  layer,  while  band  #2, 
within  the  cutoff  frequency,  will  be  excited. 
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Figure  6.4  Membrane  support  mounting  detail 
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Figure  6.6  Displacement  measurements  of  membrane  excited 

by  turbulent  boundary  layer  pressure  fluctations. 
"o",  = 32  m/sec;  "x",  Uw  = 48  m/sec. 
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Figure  6.7  Averaged  sound  pressure  level 

"o" , U = 32  m/sec.  "x",  U = 48  m/s ec. 
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Figure  7.1  Radiated  sound  power  level  estimates. 

"o",  based  on  wall  pressure  data; 

"x",  based  on  displacement  measurements; 

based  on  sound  pressure  level  measurements 


Displacement  level  in  10  Hz  band  width 


Figure  8.Lc  Displacement  level  in  10  Hz  band  width 
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FIGURE  8.2  Ratio  Between  Pitching  and  Heaving 
Resonant  Frequencies  for  a Bar 
Resting  on  Two  Strings. 
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FIGURE  8.3  String  and  Block  Model  Excited  Below 
Foundamental  Resonant  Frequency  of  a 
2. .•  x Membrane  Unit. 


